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In a mesoscopic system, under zero bias voltage, a finite charge is transferred by quantum adi¬ 
abatic pumping by adiabatically and periodically changing two or more control parameters. 
We obtained expressions for the pumped charge for a ring of three quantum dots (QDs) by 
choosing the magnetic flux penetrating the ring as one of the control parameters. We found 
that the pumped charge shows a steplike behavior with respect to the variance of the flux. 
The value of the step heights is not universal but depends on the trajectory of the control 
parameters. We discuss the physical origin of this behavior on the basis of the Fano resonant 
condition of the ring. 


1. Introduction 

Recently, quantum adiabatic pumping has attracted much attention . 1 ^ 6 Quantum adiabatic 
pumping induces a finite current in a mesoscopic system at zero bias voltage by adiabatically 
and periodically changing two or more control parameters. Classical pumping based on the 
Coulomb charging effect such as that in single-electron transistors 7 - 8 or turnstile devices 9 
does not require phase coherence. In contrast, quantum pumping is fundamentally different 
from classical pumping. Periodic deformation of two or more parts of the potential induces 
phase-coherent redistribution of the electron charges in an open quantum system. During 
this redistribution, electrons can be coherently pumped from one lead to the other leads. In 
spatially periodic systems, Thouless showed quantized charge transport induced by adiabatic 
and periodic changes in the potential . 10 Later, a formulation using a scattering matrix ap¬ 
peared, 11 - 12 which is called Brouwer’s formula. This formula is most conveniently applied to 
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a quantum adiabatic pumping in noninteracting systems. Relative modulating phase of the 
control parameters (this phase is not a scattering phase, which will appear in the following) 
in Brouwer’s formula determines the magnitude and the sign of the pumped charge. An ex¬ 
perimental demonstration of Brouwer’s formula has been reported. 13 In this experiment, two 
gate voltages controlled the periodical deformations of the shape of the quantum dot (QD), 
and then the pumped current was observed, where the amplitude of the current changed with 
the relative phase of the two gate voltages, as the theory predicted. However, the result is still 
open to argument since the pumped current can also be explained by the rectification effect 
of the displacement currents generated by the time-dependent gate voltages. 14-18 

In addition to the experimental studies, there have been several detailed theoretical studies 
on Brouwer’s formula. For example, the maximum value of the pumped charge per cycle 
becomes exactly an elementary charge 19 by appropriately choosing two potentials as control 
parameters in one QD system. The effect of the resonance on the quantum adiabatic pumping 
has been analyzed in a double-barrier quantum well 20 and a QD in a turnstile geometry. 21 The 
effect of dephasing has also been studied. 22 - 23 In two-terminal systems, the scattering matrix is 
given by a 2 x 2 unitary matrix, and there are four independent real parameters in this matrix. 
Avron et al. clarified the roles of each parameter in the transport process. 24 Additionally, the 
inverse process, namely, adiabatic quantum motors driven by applying finite bias, has been 
analyzed with scattering matrix formalisms. 25 

As discussed above, there have been many studies on quantum adiabatic pumping with 
modulating potentials. However, there have been no studies on choosing a scattering phase 
as one of the control parameters. Avron et al. considered the role of the phase in quantum 
transport in a general framework. 24 However, there seem to be no studies based on an ex¬ 
plicit model system. On the other hand, the quantum mechanical phase plays an important 
role in the field of quantum transport. The Josephson effect 26 is one of the examples, where 
the current occurs owing to the difference in the phase of macroscopic wave functions. The 
electron phase is a quantum mechanical value, and electrons obtain the phase by transport 
in the scattering region, e.g., QDs. Choosing the phase as one of the control parameters is 
physically interesting because the phase can be related to the local bias voltage between QDs 
(e.g., Faraday’s law of electromagnetic induction). 27 The effect of the time-dependent vector 
potential on the electrons in a metallic system has been explored. 28-30 Also, quantum adi¬ 
abatic pumping using the ac Josephson effect has been proposed. 31 To investigate how the 
periodicity and quantum property of the phase appear in the transport process, we analyze 
quantum adiabatic pumping by choosing a tunnel phase as one of the control parameters. 
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We treat a three-QD ring, where the modulation of the tunnel phase is equivalent to the 
modulation of the Aharonov-Bohm (AB) flux penetrating through the ring. The effect of the 
nonadiabatic modulation of the AB flux on the transport of the ring has been studied. 32-34 
Moreover, quantum adiabatic pumping with a single control parameter under a finite AB flux 
has been examined. 35 If we choose one QD energy level and the flux penetrating through the 
ring as control parameters, the pumped charge becomes a sinusoidal form or a steplike form as 
a function of the variance of the phase, depending on the values of the system parameters (the 
energy levels of QDs, the tunnel strengths). For the parameters realizing steplike behavior, 
there is no upper bound of the pumped charge. 

The structure of this paper is as follows. In Sect. 2, we explain theoretical model and 
formal solutions. In Sect. 3, we numerically analyze quantum adiabatic pumping in a three- 
QD ring, choosing the energy level of one QD and the flux penetrating through the ring as 
control parameters. In Sect. 4, we analytically discuss the origin of the singular behavior of 
the kernel under the condition of weak couplings to the leads. Finally, in Sect. 5, we conclude 
the paper. In the Appendix, we provide the derivation of the kernel of the three-QD ring. 

2. Model 

2.1 Hamiltonian 

In this section, we explain the model considered in this study, where a QD system is 
coupled to two leads, as shown in Fig. 1. The left and right leads couple to QD1 and QD2, 
respectively. We considered only a single level in each QD. We disregarded the spin degree 
of freedom of electrons and ignored inter-QD Coulomb interactions. Because we treated only 
a single level in each QD and did not consider the spin degree of freedom, then, according to 
the Pauli exclusion principle, the occupied electron number is 0 or 1; therefore, we did not 
need to consider intra-QD interactions. 

A A A A A. 

The Hamiltonian H = Mjds + Meads + M consists of the QD system part Hq D s , the non¬ 
interacting lead part M ea ds, and M T representing tunnel couplings between leads and QDs: 

Hqds - 
Meads — 

m t = 

Here, i and j are the indices of the QDs. d\ Id ,I is a creation (annihilation) operator of a 
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Fig. 1. (Color online) Schematic of the three-QD ring. 


*2 



Fig. 2. (Color online) Trajectory of two control parameters. 


localized electron in the /th QD. e, is the energy level of the /th QD. The tunnel couplings 
between the /th and /th QDs, t t j (= t*X have a Peierls phase 0,- ; - as t t j = \t i} \ . C' ak \C a iA is 

a creation (annihilation) operator of an electron of energy s a k with wave number k in lead a. 
V frk is the tunnel coupling amplitude. We examine the model at zero bias voltage and under a 
zero-temperature condition. 

2.2 Method 

Let X x and X 2 be two independent control parameters. We change the two parameters X\ 
and X 2 on the rectangular trajectory shown in Fig. 2 along the arrow. The expression for the 
pumped charge per cycle was introduced by Biittiker et al. n and formulated by Brouwer, 12 
which is given by an oriented surface integral of the kernel n over the area depicted in Fig. 2 
in the X\-X 2 phase space, 


Q_ 

-e 


= q = - 



dX ] dX 2 U(X ] ,X 2 ), 


(4) 
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where q indicates the number of pumped electrons from lead R to lead L per cycle and e is 
an elementary charge. In the following, we call q the pumped charge. II (X!, X 2 ) is the kernel 
expressed by the scattering matrix S, which is a function of two control parameters, 


U(X u X 2 ) = -Im 
n 


dS* LL dS u dS* rP dS 


+ 


LR 1 


1 LR 


dXi dX 2 ' dX { dX 2 1 ' 

The components of the scattering matrix are given by retarded Green’s functions through the 
Fisher-Lee relation, 36-38 

S ll = l~iT L G r n (s F ), (6) 

Slr = (7) 

T a (a = L or R) is the line width defined by r a = 2 np a | VJ 2 with the wide-band limit (ignoring 
the energy dependence of the line width), where p a is the density of states of lead a. The 
retarded Green’s function of the QD system, G r nm (s), is defined by the Fourier transformation 
of 


G r nm (t, 0 = - l -9(t-t')({d n (t),dl(t ')}), 


( 8 ) 


where n, m = 1,2,3 and 6 (l) is the step function. The operators are in the Heisenberg picture 
[e.g., d n (!) = e lHt/n d n e~ lHt/fl ] and (• • •) denotes the quantum mechanical and statistical average 
at zero temperature. Since we consider zero bias voltage and zero temperature condition, the 
energy of the incident electrons, s, is set to the Fermi energy s F . Hence, the kernel can be 
obtained with the retarded Green’s function: 


mxi,x 2 ) 


-Iin <! F 
n 


dG\\ (e F ) dG’ u (s F ) 
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Using the equation of motion method, the Fourier transform of the retarded Green’s function 
of the three-QD ring is given by 
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We derive the kernel Eq. (9) choosing the energy of QD 3, e 3 , and the AB phase 0 A b from the 
magnetic flux penetrating through the ring as control parameters. Details of the calculation 
are given in the Appendix. To make the notation simpler, we introduce the dimensionless 


parameters X\ 


r L -,*2- rR , ^ - y - yj—, — 

we have introduced a positive constant y to normalize the energy of QD 3 £3 (the choice of y 
does not affect the result). We also introduce the normalized kernel n(.r 3 , 0ab) = yn(e 3 , c/Ab), 
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and the pumped charge is obtained by 


r>^3max f" 

= - I dx 3 I 

^'^3min ^ 0 / 


*0ABmax 

0ABmin 


d(p AB nU'3, <p AB )- 


(ID 


We can see that the kernel Eq. (A-6) is a periodical function of (p AB . Regarding the normalized 
kernel ft (x 3 , <p AB ) as a function of (x x , x 2 , 1 ^23 UI ^31 U <P\ b)> we found the following symmetry: 


ft(-*T 


ja.jc2.h23l.h31l 


= ft(*3,' 


JC2.JC1.h31l.h23l 


( 12 ) 


which is related to the symmetry of time and space reversal. Moreover, regarding the kernel 
ft (jc 3 , (Pab) as a function of (x \, x 2 , * 3 , <pxs>), there is another symmetry, 


ft (-X 3 , 0 AB )| V1 = -ft(-JC3,^'-0 A B)|_ Jclf-Jft . (13) 

which is related to the electron-hole symmetry. 


3. Numerical Results 

In this section, we present the numerical results of the kernel and pumped charge Q. In 
Figs. 3 and 4, we show the contour plots of the kernel ftu 3 ,0 AB ) and the corresponding 
pumped charge q as a function of the upper bound of the phase 0 ABmax . The integration region 
of the AB phase is -n < (p AB < <fi ABmax . We set the normalized tunnel couplings to |.V| 2 | = 
ta 3 | = tail = 1. Figure 3 corresponds to normalized energies of QDs 1 and 2 being the 
same ( x\ - x 2 = 0.1), and Fig. 4 corresponds to normalized energies of QDs 1 and 2 being 
opposite (a'i = -x 2 = 0.1). In Figs. 3(c) and 4(c), the solid lines represent the integration 
region -3 < x 3 < 3 of the kernel and the dashed lines represent the integration region -6 < 
x 2 < 6. Positive and negative peaks appear (almost) periodically in the direction of the phase 
in the contour plot of the kernel. When the signs of the energies of QDs 1 and 2 are the 
same, the peak heights are smaller than the depth of the negative peaks (Fig. 3), and the 
pumped charges become increasingly negative with 0 ABm ax on average. When the signs of 
the energies of QDs 1 and 2 are opposite, the peak heights and depths are the same (Fig. 4). 
Then the pumped charge per cycle is periodical as a function of 0 ABmax and a sinusoidal form. 
We can understand this behavior for xi = -x 2 = 0.1 by considering the two symmetries, 
Eqs. (12) and(13), ft (* 3 , 0a B )| V 1jA ., = ft(x 3 , -(/> A b)\ X2Xi = -fl(-x 3 ,(p AB + ^)\_ X2 _ Xi (note that 
1 12 1 = IS 23 I = tail = 1)- In Figs. 3(d) and 4(d), the solid lines are the results for the integration 
region 0 < x 2 < 3 of the kernel, and the dashed lines are for the integration region 0 < a 3 <6. 
When we choose the integration region to pick up only negative peaks, the corresponding 
pumped charge shows a steplike form as a function of 0 ABma x- There is no upper (or lower) 
bound of the pumped charge per cycle, except for the special situation explained above. 
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Fig. 3. (Color online) (a) Contour plots of normalized kernel ft (* 3 , 0ab) and (b) transmission probability of 
the three-QD ring (the horizontal axis is the normalized energy of QD 3 and the vertical axis is the phase 
<Pab in unit of 2n.) The color bar indicates the value of the kernel and the probability. We set normalized tunnel 
couplings |si 2 | = [S 23 I = 1^31 i = 1. The normalized energies of QDs 1 and 2 are the same ( x\ = x 2 = 0.1). 
Arrows in (a) indicate the trajectory of the parameters, (c) and (d) Pumped charges (the horizontal axis is the 
upper bound of the phase Tahiti ax and the vertical axis is the pumped charge q)\ in (c) solid lines represent the 
integration region -3 < x 3 < 3 of the kernel, and dashed lines represent the integration region -6 < x 3 < 6 ; in 
(d) solid lines represent the integration region 0 < x 3 < 3 of the kernel, and dashed lines represent the integration 
region 0 < x 3 < 6 . 


We found very different behavior of the kernel at a special choice of the parameters. In 
the following, we restrict ourselves to the symmetric situation |s 12 | = taal = tail = Z. When 
X] = x 2 = Z, there are a series of isolated dips at jc 3 = Z and (p A b = 2/m with an integer n. 
Figures 5(a) and 5(b) show the contour plot of the kernel fl (x 3 , </> A is) and the transmission 
probability under the condition Z = 1, respectively. Figure 5(c) shows the pumped charge. 
The solid line represents the integration region 0.7 < x 3 < 1.3, and the dashed line represents 
the integration region 0 < xj < 2. Because of isolated structures in the kernel, the pumped 
charge per cycle becomes steplike as a function of the upper bound of the phase ^ABmax- This 
behavior means that the pumped charge is quantized and is robust against the uncertainty of 
the phase. Using the symmetry Eq. (13), we have a series of isolated peaks at cp AB = (2 n + 1 )n 
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Fig. 4. (Color online) Similar plot to Fig. 3 except for the normalized energies of QDs 1 and 2 being the 
opposite sign {x\ - —X 2 — 0.1). 


(, n : integer) and jc 3 = -Z when xi = X 2 = -Z. 

4. Discussion 

The origin of these peaks/dips of the kernel can be understood as the resonance behavior 
of the scattering matrix. The transmission probability of this three-QD ring is 

T= -^1*3 \s l2 \- |s23*ul e-**»\ 2 , (14) 

|a 3 |- 

where the factor A 3 is defined in the Appendix by Eq. (A-7). The kernel also depends on this 
factor ft oc |A 3 r 4 . Therefore, both the transmission probability and the kernel are enhanced 
when |A 3 | is strongly suppressed (quasi-resonant condition). We rewrite this factor as 

A 3 = [x\X2 - \s\2? — — + 21^12^23^311 COS 0AB ~ (l^l^l + l^ll 2 ^) 

+ 2 [(^1 + *2)X3 - IS 23 I 2 - |5 3 i| 2 ] • (15) 

On the basis of the numerical results shown in the previous section, we consider the condition 
of weakly coupled triple QDs to the leads, namely, Is^U^U^il » 1. Let us discuss the 
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Fig. 5. (Color online) (a) Contour plot of the kernel ft (* 3 , </>,\h) and (b) transmission probability T (the hori¬ 
zontal axis is the normalized energy of QD 3 X 3 and the vertical axis is the phase 0 ab)- (c) Pumped charge (the 
horizontal axis is the upper bound of the phase 0ABmax and the vertical axis is the pumped charge q). The solid 
line represents the integration region of 0.7 < X 3 < 1.3, and the dashed line represents that for 0 < x 3 < 2. We 
set the normalized tunnel coupling and the energies of QDs 1 and 2 to = l.s' 23 1 = [S 31 [ - X\ - xo - 1. 


n 



<Pab 


Fig. 6. (Color online) Kernel ft as a function of <f >ab at [S 12 I = 3 and [3 — 1. 


origin of the sharp isolated dip found for all positive xi,x 2 , and * 3 . In order to obtain isolated 
dips, A 3 should be strongly suppressed at <p AB = 2nn (When x\ , x 2 , and x 3 are all negative, 
resonance occurs when (p AB = (In + 1 )n). Then we require that the third term in Eq. (15), 
which is independent of jc 3 , becomes 2 \s i2 S 23 S 3 i\ and we have 


*2 = 


1^311 


^(2fti 2 523^3l|-|^23| 2 ^l), 


By substituting this relation into the first term, the coefficient of jc 3 becomes 

,2 


7 1 

*1*2 - ftd - T = 


$23 


■S31 


*1 ~ IS12I 


Therefore, the real part of A 3 becomes very small at x\ 


^12^31 


S23 
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and x 2 


S 12 S 23 


S31 
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condition is achieved at the value of jc 3 where the imaginary part of Eq. (15) becomes zero: 


(*1 + x 2 )x 3 
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^23 5 31 


512 


(17) 


Figure 5(b) shows the transmission probability as a function of the two control parameters. 
Clearly, near the condition of a sharp dip of the kernel, the transmission probability is strongly 
enhanced, but becomes zero exactly at the condition of the sharp dip of the kernel (in Fig. 5, 
v 3 = 1,0 ab = 0,2tt). This sharp dip of the conductance near the resonant conditions can 
be understood as Fano resonance , where the continuum spectra of the leads and the discrete 
energy level of the three-QD ring interfere. 

Then we introduce the offset variable a to study the behavior of the kernel near the res- 


* 12*31 

*23 


and x 2 


* 12*23 

*31 


. We estimate the 


onant dip, defined as a = x 3 - yy 1 and fix x\ 
transmission probability and the kernel for ta 2 |, |.s' 231 , |.v 3 11 » 1 , with neglecting a small term, 
namely, the factor -1/4 in the first bracket of Eq. (15) times a. The transmission probability 
reads 
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where F (/)AI , = 2\s l2 \ sin -y, and hence Fano resonance is expected near a 
Similarly, the kernel is 
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, and we have neglected the terms of 


the first order of a in the numerator since it means no contribution to the pumped charge when 
integrated in the range -a max < a < a max . In the following discussions, we restrict ourselves 
to the symmetric situation, l^l = tail = P\s\ 2 \, then the kernel becomes simpler, 

\s 12 \ (3 + F L) (1 - F L) cos <^b - (f f sin 2 
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and the kernel has simple dips at a = 0 , whose width is relatively broad as £2 
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0.32/3 2 . Ata = 0, 


fi(0, 0 ab) 


l^l2l(l- F L) COS ^B 


( 22 ) 


(i + F 2 ) 

\4 ^ r <j> ab) 

This kernel as a function of is shown in Fig. 6. Clearly, the kernel has a sharp dip at 
0 ab = 0 but changes its sign at F l/)A]i = ±1/2, namely, cp A b = The width of the dips 

becomes very narrow as IS 12 I '• 

Using the obtained kernel, we can evaluate the pumped charge for the areas, [ci min , a max ] = 

[-00, 00] and [- 0 ab max? 0 AB max]? 


q (^ab max) 


ki2l 


d(p AB 


0ABir 


(z - (2|s 12 | sin £f -) 2 ) cos cp AB 
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2 0AB 


a + (2|s 12 |sin^) 2 ) 2 


T + (2|^ 12 | sin %s.) 2 


(23) 


We then obtain analytical expressions for two extreme settings: 

• ^ABmax = 1 /(2 |.V| 2 |), where we take only the negative part of the dip, then 
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"2|s 12 I 
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(z + (2|S 12 | Sin %£) 2 ) 2 2 + (2|5 12 | sin ^) 2 


+ o 
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1^121" 


(24) 


for |5 , i 2 | —» oo. Therefore, the pumped charge per cycle is unity, which seems to corre¬ 
spond to the preceding result studied in a one-QD system with two potentials chosen as 
control parameters. 19 


(p ABmax = n - where we take both the negative and positive contributions around the dip, 


q(n) = 


n 


2\s n \ 


(16[5 12 | 2 + 1) 3/2 - 24|s 12 | 2 - 1 
(16|5 12 | 2 + l) 3 ' 2 


1 


Vl6|^l 2 | 2 + 1 


T7~\ ( for l^l °°)- (25) 

4|s 12 | 

Therefore, we do not expect the pumped charge when choosing a large area in the phase 
space since the positive and negative kernels cancel each other. 


When [a min , a max ] = [~b,b\ for finite b > 0 and [-0ABmax, 0ABmax], we could not obtain 
analytical results. Figure 7 shows a numerical estimation of the pumped charge as a function 
of 0ABmax for various values of b of 1,10,100 and |si 2 | = 10. Except for a very small (p ABmia , 
the pumped charge is nearly constant, which enables the clear steplike behavior described in 
the previous section. However, the step height is not universal and depends on the range of 
the integration [—b, b\. 
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0ABmax 

2 77 - 


Fig. 7. (Color online) Pumped charge as a function of the varying range [— 0ABmax, 0 a b max] and \-b, b] for 
b - 1 (dashed), 10 (dot-dashed), and 100 (solid) for [S12I = 10 (top) and ft — 1. 


5. Conclusions 

We obtained explicit expressions for the pumped charge for a ring of three noninteract¬ 
ing quantum dots (QDs) using Brouwer’s formula. We chose the energy of one QD and the 
Aharonov-Bohm phase determined by the flux penetrating through the ring (£ 3 ,0 A b) as con¬ 
trol parameters. We found that the pumped charge per cycle shows quasisinusoidal or steplike 
behavior depending on the energies of QDs 1 and 2, the tunnel couplings, and the variable 
range of the parameters. The step height is not universal and depends on the trajectory of the 
parameters. Explicit analytical expressions are obtained for the situation that the three QDs 
are weakly coupled with the leads, where the steplike behavior is related to the Fano reso¬ 
nance. For realizing a current standard, our results have an advantage since we can obtain 
large charges per cycle by changing the phase indefinitely. Moreover, from the viewpoint of 
the stability of the current standard against the fluctuations of the flux in the experiments, the 
steplike behavior found in the three-QD ring may enable the precise control of the pumped 
charge. 
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Appendix: Kernel and pumped charge of three-QD ring 

We consider the kernel of three-QD ring n(Aj, X 2 ). By definition, the retarded Green’s 
function has the property G r (G r ) 1 = 1. By differentiating this identity with respect to X, we 
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obtain 


dG r 

~dX 


G'fyG', 


(A-l) 


where we define the matrix f x = d (G r )“ 1 /dX. The kernel requires the (1,1) and (1,2) com¬ 
ponents of Eq. (A-l). 

When we choose X = e 3 , 
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then 

= ( ?12 ^21 “ t2iG r l2 ), j = i{tuG r n G r 2 2 - ?2i (G r n y) • (A-5) 

Although we have chosen the tunneling phase cpn as one of the control parameters, we 
found that the kernel depends only on the total sum of the tunneling phases: cp = (pn+Qn+foi ■ 
This result is invariant even if we choose other tunneling phases cp 2 3 or cp 3l as the control 
parameter. This is related to the electron coherence in the three-QD ring and gauge invariance. 
If one of the tunneling amplitudes t 2 3 or / 3] is zero, the phase coherence is not maintained 
along the ring, and the kernel may explicitly depend on <p\ 2 . We can divide the phase (p into 
two parts, the AB phase from the magnetic flux penetrating through the ring c/Ais {(Par = 27 rA-. 
©0 = -) and the phase independent of the flux cp*\ (p = 0 AB + cp*. Assuming that the QDs are 
made of 5-orbitals, tunnel coupling provides lower energy for a symmetric orbital without 
nodes, and the tunnel amplitude should take a negative value t\ 2 = - \tn\- This corresponds 
to e“ f ’* = (-1) 3 = -1, i.e., (p* = n. The phase that we can change freely is cp ab from the flux 
penetrating through the ring and we regard <p AB as a control parameter; thus the final form of 
the kernel is 


fl fe, <^ab) 
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where A 3 is the determinant of the retarded Green’s function given in Eq. (10) and given by 


+ T (*2 + r 1*3 + 2 1 ^ 12 ^ 23^311 COS0AB _ l^l” X 1 + X / I ^ 3 11 *2 + » _ 1 ^ 121 “ *3- 
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